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Total  Positivity:  A  Review 
by 

Jee  Soo  Kim  and  Frank  Proschan 

'  ABSTRACT 

This  paper  is  an  invited  entry  for  the  Encyclopedia  of  Statistical 
Sciences,  edited  by  N.  L.  Johnson  and  S.  Kot;  and  published  by  John  Wiley 

»  .  ^  *  i 

and  Sons.  The  main  objective  is  to  review  the  concepts  of  total  positivity 
which  plays  an  important  role  in  various  domains  of  mathematics  and 
statistics.  This  article  describes  the  power  and  scope  of  total  posi¬ 
tivity,  and  samples  the  great  variety  of  fields  of  its  applications.  1 " 

1.  Introduction. 

The  theory  of  total  positivity  has  been  extensively  applied  in  several 
domains  of  mathematics,  statistics,  economics,  and  mechanics.  In  statistic 
totally  positive  functions  are  fundamental  in  permitting  characterisations 
of  best  statistical  procedures  for  decision  problems.  The  scope  and  power 
of  this  concept  extend  to  ascertaining  optimal  policy  for  inventory  and 
system  supply  problems,  to  clarifying  the  structure  of  stochastic  processes 
with  continuous  path  functions,  to  evaluating  the  reliability  of  coherent 
systems,  and  to  understanding  notions  of  statistical  dependency. 

In  recent  years  Samuel  Karlin  has  made  brilliant  contributions  in 
developing  the  intrinsic  relevance  and  significance  of  the  concept  of  total 
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positivity  to  probability  and  to  statistical  theory.  In  1968,  Karlin  wrote 
a  classical  book  devoted  to  this  vast  subject.  This  remarkable  book  presents 
a  comprehensive,  detailed  treatment  of  the  analytic  structure  of  total ly 
positive  functions,  and  conveys  the  breadth  of  the  great  variety  of  fields  of 
it  applications.  This  book,  together  with  Karlin's  other  fundamental  papers, 
inspired  many  new  developments  and  discoveries  in  many  areas  of  statistical 
applications.  Frydman  and  Singer  r8'  obtained  a  complete  solution  to  the 
embedding  problem  for  the  class  of  continuous-time  Markov  chains:  The  class 
of  transition  matrices  for  the  finite  state  time- inhomogeneous  birth  and 
death  processes  coincides  with  the  class  of  non-singular  totally  positive 
stochastic  matrices.  Keilson  and  Kester  r21]  employed  total  positivity  to 
characterize  a  class  of  stochastically  monotone  Markov  chains  which  has  the 
property  that  the  expectation  of  unimodal  functions  of  the  chain  is  itself 
unimodal  in  the  initial  state.  To  help  unify  the  area  of  stochastic  com¬ 
parisons  Hollander,  Rroschan  and  Sethuraman  ‘9’  introduced  the  concept  of 
functions  decreasing  in  transposition  (DT) .  In  the  bivariate  case,  a  func- 


tion 

f(\r 
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Xj,  x,j  is  said  to  have 

the  DT 

property 

if 

(a) 
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;  x2,  Xj 

!  :  i .  e .  , 

transposing  from 

the  natural  order  ( x ^ ,  x,)  to  (x-,,  Xj)  decreases  the  value  of  the  function. 
In  their  paper,  total  positivity  is  essential  in  showing  that  i'',lii  =  r'.  the 
probability  of  rank  order  is  a  DT  function. 

Karlin  and  Rinott  r  18"!,  i  19]  extended  the  theory  to  multivariate 
cases.  Multivariate  total  positivity  properties  are  instrumental  in  IS’ 
and  [19]  for  the  results  which  are  applied  to  obtain  positive  dependence 
of  random  vector  components  and  related  probability  inequalities. 


fl 
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For  an  excellent  global  view  ot"  the  theory,  as  well  as  for  pertinent 
references,  the  reader  may  consult  Karlin  T13]. 


2.  Definition  and  Basic  Properties. 

2.1.  Definition  of  Totally  Positive  Function. 


A  function  f(x,y)  of  two  real  variables  ranging  over  linearly  ordered  one¬ 
dimensional  sets  X  and  Y,  respectively,  is  said  to  be  totally-  positive  of  order 

k  (TP.)  if  for  all  x.  <  x_  <  ...  <  x  ,  v  <  v  <  ...  <  v  (x.  in  X:  y.  in  Y), 

—  1  2  m  •  1  *  2  -mi  -i 

and  all  1  <  m  <  k, 
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Typically,  X  and  Y  are  either  intervals  of  the  real  line  or  a  countable  «et  of 
discrete  values  on  the  real  line,  such  as  the  set  of  all  integers  or  the  set 
of  nonnegative  integers.  Mien  X  or  Y  is  a  set  of  integers,  the  term  "sequence" 
rather  than  "function'  is  used.  If  f(x,y)  is  TP.  for  all  positive  integers 
k  =  l,  2 . then  f(x,y)  is  said  to  he  totally-  positive  of  order  written 


TP  or  TP. 

oo 

A  related,  weaker  property  is  that  of  sign 
is  sign  regular  of  order  K  (SR^l  if  for  every  x 

<  v  ,  and  1  <  m  <  k,  the  sign  f 
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Many  well  known  families  of  density  functions  (both  continuous  and  discrete! 

are  totally  positive.  It  should  be  noted  that  TP?  is  the  order  of  TP-ness 

which  has  found  greatest  application.  In  the  context  of  statistics,  the  TP. 

property  is  referred  to  as  the  monotone  likelihood  ratio  property.  Higher 

order  TP-ness  is  hardly  used  in  application  except  for  the  occasional  use  of 

TP.. 

a 

Some  examples  of  functions  that  possess  the  TP  property  are: 

( i T  f(x,y)  =  exy  is  TP  in  x,ye(-®,°°),  so  that  f(x,y)  =  x^  is  TP  in 
and  yeC-00,00) . 

(ii)  f(k,t)  =  e  At[(At)k/k!]  is  TP  in  te(0,»>)  and  k  {0,  1,  2,  ...'. 


(iii)  f (x,y) 


1  it'asxsysb 
0  if  a  <  y  <  x  <  b. 


2.2.  PF^  as  Special  Case  of  Interest. 

The  concepts  of  TP^  and  TP.  densities  are  familiar  ones,  livery  density  is 
TP^j  while  the  TP?  densities  are  those  having  a  monotone  likelihood  ratio. 

A  further  important  specialization  occurs  if  a  TP^  function  may  be  written 
as  a  function  f(x,v)  =  t'(x-v)  of  the  difference  of  x  and  v,  where  x  and  y 
traverse  the  entire  real  line;  f(u)  is  then  said  to  be  a  Polya  frequency  function 
of  order  k(PF^,)  .  Note  that  a  P<51ya  frequency  function  is  not  necessarily  a 

/OO 

probability  frequency  function  in  that  j  ^ t(uldu  need  not  be  1  nor  even  finite. 

The  class  of  PF?  functions  is  particularly  important  and  has  rich  appli¬ 
cations  to  decision  theory  [10],  [11],  [12],  [18],  reliability  theory  '5],  and 
the  stochastic  theory  of  inventory  control  models  LI],  Tib]. 

Every  PF7  function  is  of  the  form  e  where  ( x]  is  convex.  On  the 

other  hand,  there  exists  no  such  simple  representation  for  PF^,  k  e  3.  Prob¬ 
ability  densities  which  are  PF  ,  abound.  For  other  properties  and  examples 
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of  PF0  densities,  see  the  entry  "Pdlya  Type  2  Frequency  Distributions." 

Probability  densities  which  decrease  to  zero  at  an  algebraic  rate  in 
the  tails  are  KOt  PFT,.  For  example,  (i)  Weibull  with  shape  parameter 
<  1:  f(x)  =  aX(Xx)a  ^exp[ - ( Xx)a] ,  x  >  0,  X  >  0,  0  <  a  <  1,  and  (ii'j 
Cauchy:  f  ( x  1  =  1/l’tt  (l+x“)  ] ,  <  x  <  <*>  are  not  PF\. 

Intriguing  results  in  the  structure  theory  of  PF^  functions  can  be  found 

in  Karlin  and  Proschan  [16],  Karlin,  Proschan,  and  Barlow  [17],  and  Barlow 
and  Marshall  [2]. 

2.3.  Variation  Diminishing  Property. 

An  important  feature  of  totally  positive  functions  of  finite  or  infinite 
order  is  their  variation  diminishing  property:  If  ffx.y)  is  TPj.  and  giyl 
changes  sign  at  most  j  <  k-1  times,  then  h(x)  -  f  ffx.y)  g(y)  dy  changes  sign  at 

i 

most  j  times;  moreover,  if  h(x)  actually  changes  sign  j  times,  then  it  must  i 

i 

change  sign  in  the  same  order  as  gfy)  .  It  is  this  distinctive  property  which 
makes  TP  so  useful.  The  variation  diminishing  property  is  essentially  equiva- 

I 

lent  to  the  determinantal  inequalities  (1).  Greater  generality  in  stating 
this  property  is  possible.  The  interested  reader  is  referred  to  Chapter  1, 

Karlin  [13].  A  more  direct  approach  to  the  theory  is  taken  by  Brown  et  al.  r-  5 ' , 
giving  appropriate  definitions  and  criteria  for  checking  directly  whether  a 

*  i 

family  of  densities  possesses  variation  diminishing  property. 


| 

I 

2.4.  Composition  and  Preservation  Properties. 

Many  of  the  structural  properties  of  TP^  functions  are  deducible  from 
the  following  basic  identity  which  is  an  indispensible  tool  in  the  study  of 
j  total  positivity. 

|  Basic  Composition  Formula.  Let  h(x,t)  =  {  f(x,y)  gly.t]  do(y]  converge  abso- 
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lutely,  where  do(y)  is  a  sigma-finite  measure.  Then 


x  |  > x  2  >  • 

"Xnl 

xl’x2; 

t1,t2,  . 

1  =  /  /  fj 

1  ' 

-v-v 

> 

>Vy2’  ' 

•g 

j 

*1 ’^2  ’  ' 

J 

da(v^) 

A  direct  consequence  of  the  composition  formula  is:  If  f(x,v)  is  TPffl 

and  g(y,t)  is  TP^,  then  h(x,t)  =  (  f(x,v)  g(y,t)  defy)  (the  convolution  of  f 

and  g)  is  TP  .  ,  . .  In  manv  statistical  applications  this  consequence  is 

exploited  principally  in  the  case  when  f  and  g  are  Polya  frequency  densities. 

That  is,  if  f(x)  is  PF;m  and  g(x)  is  PFn>  then  h(x)  =  /  f(x-t)  g(t)dt  is 

PF  „->•  From  this  ive  can  obtain  a  kev  result  as  follows. 

min(m,n) 

Theorem  1 .  Let  f ^ ...  be  density  functions  of  nonnegative  random  vari¬ 
ables  with  each  f ^  a  PF^.  Then  g(n,x)  =  f^  *  f  *  fn(x)  (*  indicates 

convolution)  is  PTj.  in  the  variables  n  and  x,  where  n  ranges  over  1,  2,  ...  and 
x  traverses  the  positive  real  line. 

The  case  when  the  random  variables  are  not  restricted  to  be  nonnegative 
is  discussed  in  Karlin  and  Proschan  [161.  These  composition  and  preservation 
properties  allow  us  to  generate  other  totally  positive  functions,  thus  making 
it  easy  to  enlarge  the  TP  or  PF  classes  and  to  determine  whether  the  TP  pro¬ 
perty  holds. 

2.5.  Unimodality  and  Smoothness  Properties. 

A  function  totally  positive  or  more  generally  sign  regular  is  endowed  with 
certain  structural  properties  pertaining  to  unimodality  and  smoothening 
properties.  From  the  definition  of  PF^  can  be  derived 

i  I 


7 


ffxj-y)  -f'Cxj-y) 
f(x2-y)  -f'(x2-y) 


>  0 


f3) 


for  Xj  <  x,  and  y  arbitrary. 

In  the  event  that  f  (uq)  =  0  the  above  inequality  implies  that  f  fu)  >  0  for 

u  <  and  f  (u)  2  0  for  u  >  Uq.  This  clearly  implies  that  if  f(u)  is  PF 

then  f(u)  is  unimodal.  In  particular,  every  PF_,  density  is  a  unimodal  density. 

We  note  that  the  unimodalitv  result  is  valid  in  case  f  is  a  PF.,  sequence. 

We  now  describe  a  smoothening  property  possessed  by  the  transformation 
under  which  convexity  in  g(x)  is  carried  over  into  convexity  in  hfxl,  viz. 

h(n)  =  {  f ^ (xj g(x)dx  for  n  =  1,  2,  ...,  (4) 

where  f  (xl  is  the  n-fold  convolution  of  f.  To  make  this  notion  precise 
assume  f(x)  is  PF,  and  g(xl  is  convex.  Let  u  =  (  xf(xldx.  Note  that  for  arbi¬ 
trary  real  constants  a^  and  a ^ , 

\  {gCxl-rCty/uIx+ajJ}  f(nl(x)dx  =  h(n)  -  (.a^n+aj) .  (SI 

Since  g(x)  is  convex,  then  g(xl -[ (a^ /u)x+a^ ]  has  at  most  2  changes  of  sign 
and  if  2  changes  of  sign  actually  occur,  they  occur  in  the  order  +  -  +  as  x 
traverses  the  real  axis  from  to  +  °°.  Since  f  is  PF_,  f^fxl  is  TP.  in 
the  variables  n  and  x  by  Theorem  1.  The  variation  diminishing  property  implies 
that  h(n)  -  (a  n+a^)  will  have  at  most  2  changes  of  sign.  Moreover,  if 
h(n)  -  (a^n+a^)  has  exactly  2  changes  of  sign,  then  these  will  occur  in  the 
same  order  as  those  of  g(x)  -  [ (aQ /u)x+a^ ] ,  namely  +  -  +.  Since  a^ ,  a^  are 
arbitrary,  we  easily  infer  that  h(n)  is  a  convex  function  of  n.  Similar 
results  apply  for  concavity. 
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5.  Applications  to  Statistical  Decision  Theory. 

Historically  this  is  perhaps  the  first  area  of  statistics  benefiting 
from  the  application  of  TP  due  to  the  great  papers  of  Karlin  [10],  'll],  '12]. 
We  consider  the  problem  of  testing  a  null  hypothesis  against  its  alternative 
hypothesis,  i.e.,  a  2-action  statistical  decision  problem.  There  exist  two 
loss  functions  and  L-,  on  the  parameter  space  where  L^(9)  is  the  loss 
incurred  if  action  i  is  taken  when  9  is  the  true  parameter  value.  The  set  in 
which  L,(9)  <  (>)  L,(9)  is  the  set  in  which  action  1  (action  2)  is  preferred 
when  9  is  the  true  state  of  nature.  The  two  actions  are  indifferent  at  all 
other  points.  We  shall  assume  that  Lj(9)  -  L_,(9)  changes  sign  exactly  n  times 
at  elf  e2,  ....  en. 

Let  i  be  a  randomized  decision  procedure  which  is  the  probability  of 

taking  action  2  (accepting  the  alternative  hypothesis)  if  x  is  the  observed 

value  of  the  random  varibale  X.  Let  C  be  the  class  of  all  monotone  randomized 

n 

decision  procedures  defined  by 


1  for  x„.  < 


;i+l ’  1  =  °> 


4>(x)  =  •  for  x  =  x.,  0  $  T  j  =  1,  2 .  n  (6) 

0  elsewhere, 

where  [a]  denotes  the  greatest  integer  <  a  and  x  = 

Using  the  variation  diminishing  property,  Karlin  [11]  obtained  the  main 
results,  which  state: 

Theorem  2 .  Let  f(x,9)  be  a  strictly  TPn+j  density  and  p(9,4>)  =  J  [(1  -  <fc(x]l) 
Lj(9)  +  4>(x)  L2(9)]  f (x , 9)  du(x).  Then  for  any  randomized  decision  procedure  4 
not  in  Cn  there  exists  a  unique  4>°  such  that  f  ( 0 , d> ])  s  f(9,4>)  with  inequalities 


everywhere  except  for  9  =  9^,  9?,...,0n. 


then 


Theorem  3.  If  *  and  y  are  two  procedures  in  C  and  f  is  strictlv  TP 
j  [  5  (xl -y  (x]  }  fCx.^i  du(.x)  has  less  than  n  zeros  counting  multiplicities. 

Assume  f  ( x ,  d )  is  strictly  TP.,.  For  a  one-sided  testing  problem,  existence 
of  a  uniformly  most  powerful  level  x  test  can  be  easily  established  by  Theorem 
2  and  Theorem  3. 

More  detailed  discussions  and  other  decision  theoretic  applications  can 
be  found  in  Karlin  [10"!,  [11],  r12?  and  Karlin  and  Rubin  r20]. 

4.  Applications  in  Probability  and  Stochastic  Processes. 

Let  P(t,x,E)  be  the  transition  probability  funct  .  uf  a  homogeneous 
strong  Markov  process  whose  state  space  is  an  interv  ^n  the  real  line  and 
which  possesses  a  realization  in  which  almost  all  st  p1  paths  are  continuous. 
Karlin  and  McGregor  ! 14  !  established  the  intimate  relationship  between  the 
general  theory  of  TP  functions  and  the  theory  of  diffusion  stochastic  processes. 
Their  main  result  shows  the  transition  probability  function  Pft.x.E'  is 
totally  positive  in  variables  x  and  E.  That  is,  if  Xj  <  x,  <  ...  <  x^  and 
E^  <  E-,  <  .  .  .  <  E n ( E ^  <  Ej  denotes  that  x  <  y  for  every  x  «■:  E-  and  y  t  E.i. 
then  det  ||  P(t  ,x^  ,Ej]  ||  >  0  for  every  t  >  0  and  integer  n.  This  relation  intro¬ 
duces  the  concept  of  a  TP  set  function  f(x,E)  =  P(t,x,E)  where  t  is  fixed,  x 
ranges  over  a  subset  of  the  real  line,  and  E  is  a  member  of  a  given  sigma 
field  of  sets  on  the  line. 

If  the  state  space  of  the  process  is  countably  discrete,  then  continuity 
of  the  path  functions  means  that  in  every  transition  of  the  process  the  particle 
changes  "position",  moving  to  one  of  its  neighboring  states.  Tims,  discrete 
state  continuous  path  processes  coincide  with  the  so-called  birth-death  pro¬ 
cesses  (Karlin  and  McGregor  ["15])  which  is  a  stationary  Markov  process  whose 
transition  probability  matrix  P.j(t)  =  Pr(x(tl  =  j  |  x(0)  =  i]  is  totally 
positive  in  the  values  i  and  j  for  every  t  >  0. 


Two  concrete  illustrations  of  transition  probability  functions  that 
arise  from  suitable  diffusion  processes  are  r14’: 

li)  Let  L^(x’)  be  the  usual  Laguerre  polynomial,  normalized  so  that  l/  = 

i. nn a ")  ,  and  let  Pit)  be  the  infinite  matrix  with  elements 

P  (tl  =  |”e" Xt  L"1  f  x)  L"1  (xlxae  Xdx. 
mn  ' 0  n  m 

Then  P(t)  is  strictly  TP  for  each  fixed  t  >  0  and  i  >  -1.  This  is  an  exa-;- : 
of  a  transition  probability  matrix  for  a  birth-death  process. 

(ii)  The  Wiener  process  on  the  real  line  is  a  strong  Markov  process  with 
continuous  path  functions.  The  direct  product  of  n  copies  of  this  process  ; 
the  n-dimens  ional  Wiener  process  which  is  known  to  be  a  strong  Market  proce  - 
Therefore  the  transition  probability  function  P(t,x,fi:  =  1 ,  <  » 4 - 1  '  ^  exp 
" - ( x-y l " Mt  "!dy  is  totally  positive  for  t  >  0. 

5.  Applications  in  Inventory  Problem. 

Suppose  that  the  probability  density  ffx)  of  demand  for  each  period  is 
PF_.  The  policy  followed  is  to  maintain  the  stock  size  at  a  fixed  level  s 
which  will  be  suitably  chosen  so  as  to  minimize  appropriate  expected  costs, 
is  determined  by  a  fixed  capacity  restriction.  At  the  end  of  each  period  an 
order  is  placed  to  replenish  the  stock  consumed  during  that  period  so  that  a 
constant  stock  level  is  maintained  on  the  books.  Delivery  takes  place  n 
periods  later.  The  expected  cost  for  a  stationary  period  as  a  function  of 
the  lag  is 

L(n)  =  h(S-y)f(nl(y)dy  +  /”  p(y-S)f (n1  (y'ldy  i~ 

where  S  is  fixed,  h  represents  the  storage  cost  function  and  c  the  penalty 
cost  function. 


Let  h  and  o  be  convex  increasing  functions  with  h(01  =  o 1 0 1  =  0.  Then 
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vs e  may  vst'ite 

L( n)  =  f  r ( y ) f 1 n 1 ( y]dy , 

j  h(S-y)  for  0  <  y  <  S 

where  r(y)  =  ^ 

|  oly-Sl  for  S  <  y. 

Novs  rlyl  is  a  convex  function.  Hence  by  the  convexity  preserving  property  of 

i-l',  we  conclude  that  I. in)  is  a  convex  function.  I'hus ,  if  the  length  of  lac 
increases,  the  marginal  expected  loss  increases. 

Very  interesting  applications  of  total  positivity  are  found  in  system 
supply  problems.  Suppose  we  wish  to  determine  the  intitial  spare-parts  kit 
for  a  complex  system  which  provides  maximum  assurance  against  system  shutdown 
due  to  shortage  of  essential  components  during  a  period  of  length  t  under  a 
budget  for  spares  C  .  We  assume  that  a  tailed  component  is  instantly  re¬ 
placed  by  a  spare,  if  available.  Only  spares  originally  provided  may  be  used 
for  replacement,  i.e.,  no  resupply  of  spares  can  occur  during  the  period.  The 

system  contains  d^  operating  components  of  type  i,  i  =  1,  2 . k.  The 

length  of  life  of  the  operating  component  of  the  i^  type  is  assumed  to  ve 
an  independent  random  variable  with  PI-'  density  f .  .  ,  i  =  l,  2,  ....  d..  The 

K  1 J  I 

unit  cost  of  a  component  of  type  i  is  c. . 

Our  problem  is  to  find  n.,  the  number  of  spares  initiallv  provided  of  the 
k  1  k 

it  '  type,  such  that  IT  P^(ni")  is  maximized  subject  to  £  njC;  •_  c  and 
i=l  i=l 

n^  =  0,  1,  2,  ...  for  i  =  1,  2,  ...,  k,  where  P^iml  =  probability  of  experi¬ 
encing  <  m  failures  of  type  i. 

In  Black  and  Proschan  r4J,  a  detailed  discussion  of  methods  is  given  for 
computing  the  solution  when  each  £nP^(ml  is  concave  in  m,  or  equivalently, 
when  each  P^(n-m)  is  a  TP,  sequence  in  n  and  m.  To  show  P.(n-ml  is  a  TP, 
sequence  in  n  and  m,  we  note: 


1 


* 


1.  c..(n),  the  probability  of  requiring  n  replacement <  of  uncrating  component 

1,  i,  is  a  PF,  sequence  in  n  for  each  fixed  i  and  l. 

2.  c^(n),  the  probability  of  requiring  n  replacements  of  type  i  ,  is  a  PF , 

sequence  in  n  for  each  i,  since  c.(n)  =  c^  *  Cp  *  ...  *  c-^  (n). 

i 

5.  P.(n-m)  is  a  TP,  sequence  in  n  and  m  for  each  i,  since 


and 


(a) 


(b) 


F-(n)  ~  c(n-m)q(m),  where 

(  1  for  m  =  0 ,  1 ,  2 , 
q(ml  =  -j 

I  0  otherwise, 
q (m)  is  a  PF  sequence, 


(cl 


the  convolution  of  PF^,  is  PF^ 


A  procedure  for  computing  the  optimal  spare  parts  kit  in  terms  of  P.iml 

is  given  in  :41:  For  arbitrarv  r  >  0,  for  those  i  such  that  inP. ( P -InP. i  (M  • 

•  i  i 

rc^,  define  n?(r)  =  0;  for  the  remaining  i,  define  n*(rl  as  1  +  rlargest  n 

such  that  ilnPjfn+li  -  UnP^n)  >  rcP.  Compute  crn*(r)]  =  £  cm*  Ir)  . 

i  =  l 

n*  is  optimal  when  c^  is  once  of  the  values  assumed  by  c’’n*ir)'  as  r 
varies  over  (0,®). 

b.  Applications  in  Reliability  and  Life  Testing. 

A  life  distribution  F  is  said  to  have  increasing  (decreasing)  failure 
rate,  denoted  by  IFR  (DFR)  ,  if  log  fl-F(t'|]  =  log  F(t)  is  concave  (convex  on 
r0,=°)).  If  F  has  a  density  f,  then  the  failure  rate  at  time  t  is  defined  by 

r(t)  =  f(t)/F(t)  for  F(t)  <  1.  Distributions  with  monotone  failure  rate  are 
of  considerable  practical  interest  and  such  distributions  constitute  a  very 


large  class. 


The  monotonicity  properties  of  the  failure  rate  function  r i t :  are 
intimately  connected  with  the  theory  of  total  positivity.  The  statement 
that  a  distribution  T  lias  an  increasing  failure  rate  is  equivalent  to  the 
statement  that  Fi.x-y)  is  TP,  in  x  and  y,  or  Fix,  is  pp , . 

The  concept  of  IT  yields  fruitful  applications  in  shock  models.  We 
say  that  a  distribution  F  has  increasing  failure  rate  average  (  1 F'  U  A  ■  if 
ll.tl  -log  Flti'  is  increasing  in  t  •  0,  or  equivalentlv,  Jit  I  1  is  de¬ 
creasing  in  t  >  0.  An  I  FRA  distribution  provides  a  natural  description 
of  coherent  system  life  when  system  components  are  independent  1 1 R .  The 
I  FRA  distribution  also  arises  naturally  when  shocks  occur  randomly  according 
to  a  Poisson  process  with  intensity  *• .  Hie  i^1  shock  causes  a  random  amount 
X^  of  damage,  where  X^,  X,,  ...  are  independently  distributed  with  common 
distribution  F. 

A  device  fails  when  the  total  accumulated  damage  exceeds  a  specified 

capacity  or  threshold  x.  Let  H(t)  denote  the  probability  that  the  device 

survives  0 , t  A 

Then  » 

|  T  T,  e  t(\t)  / k  1  for  0it<  -v 
I  k=0 

lift)  =  {  (9i 

I 

1  1  for  t  '  0. 

V. 

Note  that  e  ( \t )  \'k !  1  represents  the  Poisson  probability  that  the  device 

experiences  exactly  k  shocks  in  10, t!,  while  P^  =  F^'ix)  represents  the 

probability  that  the  total  damage  accumulated  over  the  k  shocks  does  not 

exceed  the  threshold  x,  with  1  =  P  >  P,  >  P„  >  .... 

o  1  2 

As  key  tools  in  deriving  the  main  result  in  shock  models,  the  methods 


of  total  positivity  are  employed  and  in  particular  the  variation  diminishing 
property  of  TP  functions. 


0 


1.  has  at  most  one  si  jin 


It'  I’j/  ^  is  decreasing  in  k,  P 
change,  from  +  to  -  if  one  occurs .  Then  it  follows  from  the  variation 
diminishing  property  that  ill(t)!^  ^  is  decreasing  in  t,  i.e,,  II  is  II-RA. 

The  following  inplications  are  readily  checked: 

PF,  density  II-'R  distribution  -  I  FRA  distribution. 

For  further  discussion  and  illustrations  of  the  usefulness  of  total 
positivity  in  reliability  practices  we  refer  to  Barlow  and  Prosehan  ,5'. 

P.  Multivariate  Total  Positivity  and  its  Relationship  to  Qualitative 
Notions  of  Dependency. 

The  following  natural  generalisation  of  TP,  was  introduced  and  studied 
by  Karlin  and  Rinott  : 181. 

Definition .  Consider  a  function  f(xl  defined  on  X  =  X  X,-  .  •  Xj.  where 
each  X  is  totally  ordered,  lv’e  say  that  f(xl  is  multivariate  totally  positive 
of  order  2  or  MTP,  if 

f(x  v  vi  f  (x  a  y)  i-  f(xl  f(vi  for  every  x,  y£X,  (ID1 

where  x  v  v  =  (max(Xj,  y^l  ,  max(x,,  y,l ,  ....  max(xn>  yn i i  and 
x_  a  y  =  (min(Xj,  y^|,  min(x,,  y,t,  ....  min(xn,  yn) )  .  In  order  to  verify 
(1  it  suffices  to  show  that  f(x)  >  0  is  TP,  in  every  pair  of  variables 
where  the  remaining  variables  are  held  fixed. 

Multivariate  normal  distributions  constitute  an  important  class  of  MTP, 
probability  densities.  Let  X_  follow  the  density 

f(.x)  =  (2H1  n/“|v|  1,/“expr -'2lx-u)  'B(x-u_l 

where  E  =  B  =  ||  b.^.||?  .  ^  .  This  density  is  TP,  in  each  pair  of  arguments 
and  hence  MTP,  if  and  only  if  Lr.  <0  for  all  i  *  j. 


In  a  great  many  situations,  the  random  variables  of  interest  are  not 

independent.  To  appropriately  model  these  situations  Esary,  Proschan  and 

Wat kup  '  b '  introduced  the  concept  of  association  of  random  variables.  Random 

variables  X  j ,  ....  X  are  said  to  be  associated  if  Cov(f(X),  giX)'  i-  0 

for  all  pairs  of  increasing  functions  f  and  g. 

It  can  be  shown  IS'  that  if  X  =  (X,,  X,,  ....  X  )  has  a  ioint  MTP, 

—  1  2  n 

density,  the  1:  ?|X)v|X  j  ’  a  l:.  MXi  MX.  .  t  X  i  provided  :  and  .  are  s  imu  1 1  anovus  1  v 

monotone  increasing  (or  decreasing!.  Equivalently,  CoviilXi,  ..  i  X  i  i  _•  0.  Thus 

an  MTP,  random  vector  X  consists  of  associated  random  variables. 

It  is  well  known  that  the  union  of  independent  sets  of  associated  random 

variables  produces  an  enlarged  set  of  associated  random  variables.  Clearly 

increasing  functions  of  associated  random  variables  are  again  associated.  It 

follows  that  if  X  and  V  are  independent  random  variables  each  with  associated 

components,  then  the  components  of  Z  =  X  ♦  V  is  associated.  Thus,  in  particular, 

if  X  and  V  both  have  MTP.,  densities,  then  association  of  (I,  ,  Z,,  ....  I  1  is 

retained.  However,  Z  need  not  have  a  joint  MTP,  density. 

A  key  to  many  of  the  results  on  positive  dependence  and  probabi 1  istic 

inequalities  for  the  multinormal,  multivariate  t,  and  Kishart  distributions 

obtained  by  Karlin  and  Rinott  19;  is  the  degree  of  MTP,  property  inherent  in 

these  distributions.  Their  main  theorem  delineates  a  necessary  and  sufficient 

condition  that  the  density  of  ,X|  =  (!x^|,  |X,!,  ....  JXn  1  where 

X_  =  f X j ,  X,,  ....  Xn)  is  governed  by  N(0,  "1  be  MTP,  is  that  there  exists  a 

diagonal  matrix  D  with  diagonal  elements  *1  such  that  the  off-diagonal  elements 

of  -DE  are  all  nonnegative.  For  an  illustration  of  the  power  of  this  theorem 

consider  |x|  =  ( j  X  ^  I  ,  |X^|,  ....  jX^j)  possessing  a  joint  MTP,  density  where 

P 

X  ~  N(0,  Si.  Define  S.  =  £  XT  ,  i  =  1 ,  2 ,  . . . ,  n ,  where 

1  v  =  l  l'; 

X  =  (X  ,  \  ,  X  ),  v  =  1,  2,  ...»  p  are  i.i.d.  random  vectors 

— v  1  j  2  j  n\) 


satisfying  the  condition  of  the  theorem.  The  random  variables  Sj ,  S, . 

are  associated  and  have  the  distribution  of  the  diagonal  elements  of  a  random 

positive  definite  n> n  matrix  S  where  S  follows  the  Wishart  distribution 

K  (p,  II  with  p  degrees  of  freedom  and  parameter  :: .  It  is  established  in 
1  n 

119]  that  P  (S,  >  c,,  S0  >  Ct . Sn  2  c  )  s  n  Pf  fS  ^  >  c^j  for  any  positive 

“  i  =  l 

For  other  applications  and  ramifications  of  MTP^,  see  Karlin  and  Rinott 
[18],  [19]. 

Fahmy  et  al.  7]  exploited  the  concept  of  MTP-,  to  obtain  interesting 
results  on  assessing  the  effect  of  the  sample  on  the  posterior  distribution  in 
the  Bayesian  context. 
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